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Electromagnetic moments of quasi-stable particle 
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We deal with the problem of assigning electromagnetic moments to a quasi-stable particle (i.e., 
a particle with mass located at particle's decay threshold). In this case, an application of a small 
external electromagnetic field changes the energy in a non-analytic way, which makes it difficult to 
assign definitive moments. On the example of a spin-1/2 field with mass M« interacting with two 
fields of masses M and m, we show how a conventionally defined magnetic dipole moment diverges 
at M« = M -\- m. We then show that the conventional definition makes sense only when the values 
of the applied magnetic field B satisfy \eB\/2M* <^ |A/, — M — m\. We discuss implications of these 
results to existing studies in electroweak theory, chiral effective-field theory, and lattice QCD. 



Electromagnetic (e.m.) moments of a particle are de- 
termined through observations of the particle's behavior 
in an applied electromagnetic field. For example, the 
magnetic moment is measured by observing the spin pre- 
cession in a magnetic field. In doing so, one assumes that 
the uniform magnetic field B induces a linear response 
in the energy: 



(1) 



with fl being the magnetic moment. This method works 
perfectly well for stable particles (electron, proton), as 
well as for many unstable particles (niuon, neutron, etc.), 
which live long enough for their spin precession to be 
observed. In this letter we examine the case of a "quasi- 
stable" particle, i.e., a particle with mass M* that could 
decay into two (for simplicity) particles with masses M 
and m, such that 



A4 = AI + m. 



(2) 



It turns out that applying the magnetic field in this 
situation does not lead to a polynomial energy shift 
but to a response which is non-analytic in B, typically 
AE ~ The square-root behavior is characteristic 

for the particle- production cut. In a more general situa- 
tion, when Mh. k M + m, a polynomial expansion in B 
can be made as long as 



\B\ X [magneton] <C — A/ — ' 



(3) 



which thus becomes a condition for the magnetic moment 
to be observable. 

We do not yet know of examples in nature where the 
masses of particles would be tuned to such an extent that 
the condition Eq. ^ would be violated. For example, 
the neutron mass is less than 1 MeV above the threshold 
[Mn—Mp—nie ~ 0.8 MeV), but this number is huge when 
compared to any reasonable value of the magnetic field 
measured in units of nuclear magneton: //jv — 3 x 10~^^ 



MeV/Tesla. Nevertheless, situations where the condition 
Eq. ([3]) is violated are sometimes encountered in theoreti- 
cal studies. In the studies of the VF-boson's magnetic and 
quadrupole moments as a function of bottom- and top- 
quark masses to;, and m^, a singularity at mi, + mt — Mw 
arises from the bt (or th) loop contributions. This sin- 
gularity was reported firtstly in P, Q at a time when 
the value of rut was not known yet. In lattice Quantum 
Chromodynamics (QCD), the e.m. moments of hadron 
are computed for various values of light quark masses 
and, as calculations based on chiral perturbation theory 
show, cups and singularities arise too 0, 01- In this work 
we find that the singularities arise in the region where 
the electromagnetic moments are ill-defined, because the 
condition Eq. ([3]) is not satisfied. 

Our findings are best demonstrated on a simple toy 
model of three fields: a scalar Lp and two Dirac spinors if) 
and If', interacting via the Yukawa type of coupling: 



Ant =5 9i;ip + ilj ^ip' 



(4) 



with .9 <C 1, a small coupling constant. We denote the 
masses of ip and respectively as: to, M, and M*, 
and will later on focus on the region specified by Eq. ([2|) . 

Suppose the field If", as well as one of the other two 
fields, has an electric charge e, and couples minimally to 
electromagnetism. We look for its anomalous magnetic 
moment (a.m.m.) at leading order in the coupling g. 
Depending on whether 1^9 or -0 is charged we ought to 
consider the electromagnetic vertex corrections shown in 
Fig. [1] and obtain (unprimed: ip charged, ip neutral, or 
primed: "0 charged, ip neutral): 



2g^ 



dx 



— {r + x) x{l — x) 



(47r)2 Jo ^t^^ ^ 2;(1 — x) + (1 — x)r'- 



29' 
(47r)2 



dx 



{r + x){l-x)^ 



xji^ — x{l — a;) + (1 — x)r'^ 



(5) 



(6) 



where r = M/M^, fi = m/M^ 



2 



Figure 1: One-loop electromagnetic vertex corrections. Dou- 
ble lines, single and dotted lines denote the propagators of !/', 
^p, and if, respectively. Dots denote the Yukawa coupling and 
rectangles the minimal electromagnetic coupling. 



We have checked that for = — A/ = Mn and m = 
TOtt being respectively the mass of the nucleon and the 
pion, these expressions reproduce results of the meson 
theory (the same result also arises in chiral perturbation 
theory at next-to-leading order Q). The minus sign in 
front of M appears due to the pseudo-scalar nature of 
pion. 

At Af* = m + M (or, 1 = fx + r), the denominator 
in the integrands takes the form [x^ — (1 — a;)r]^, which 
leads to an essential singularity in these expressions for 
any positive fi and r. This can explicitly be seen, for 
instance, in Fig. [2] where plotted as a function of fi. 
If if is pseudo-scalar, AI flips the sign in these expressions, 
and the singularity is replaced by a cusp. 

The singularity is clearly unphysical, since an infi- 
nite value of the magnetic moment would correspond 
to a infinite-energy response to an external magnetic 
field. To find the correct answer we consider the self- 
energy of the ip'-field in a constant electromagnetic field, 
Ffj_,^ = d^Ai, — dyA^ = const. Calculations of this sort 
have been done before, most notably by Sommerfield and 
Schwinger 0, 0| as a technique to obtain the correction 
term of order a^^^ ~ (1/137)^ to the electron's a.m.m.. 
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Figure 2: The anomalous magnetic moment «:* of !f'-field as 
function of </3-field mass ^, at fixed value r = 0.9. The red 
(solid) curve shows the real part and the blue (dashed) curve 
the imaginary part of k. (The sign of the imaginary part is 
determined by the ie prescription.) 



To cast this technique into a modern field-theoretic lan- 
guage, we introduce the sources 6*, and j for the fields 
Tp and Lp, respectively, and write down the generating 
functional of the theory. 



Z[e,,j-A] = cxp {-g J d^z(^ 



9 
6^ 



i)} 



exp 



5j{z)S9{z)Siz 
X 9ix)Six-y;A)9{y) + . 



- / d'^xd^y 



(7) 



where 



S{x - y; A) = [z7^^ - eA^ix)^ - M*] 's^'\a 



y)- 

(8) 

is the propagator of a charged Dirac particle in the pres- 
ence of an e.m. field. We then calculate the energy shift 
induced by the iZ'-field self-energy correction in the pres- 
ence of a constant e.m. field. The dependence on the e.m. 
field comes in the form of the ^^^^F^^ structure sand- 
wiched between the free tf'-field states. When the electric 
contribution is zero, this structure simply yields the pro- 
jection of the magnetic field onto the spin direction. 

The resulting energy-shift, to leading order in g, is for 
the two cases given by: 



9 



dx (r + x) 



X In 



x{l -x)B 



xy? — x{l — x) + {1 — x)r 



AE' 



dx{r + x) 
(47r)^ Ja 



(9) 



(10) 



X In 



1 - 



{l~x)^B 



xfi^ ~ x{l — x) + {I ~ x)r'^ - 
where the following dimensionless variables are used: 



B = 



eB, 
A/2 



AE 1 
Ai7 



(11) 



with Bz the projection of the magnetic field on the spin 
direction. The quantity AE is the energy shift (in units 
of A/*) due to the a.m.m. effect. In the following we will 
discuss the unprimed contribution, the primed one can 
be obtained analogously. 

The e.m. field is assumed to be small in comparison 
with the mass-scale of particles, and therefore some terms 
which are higher-order in B^ can be neglected. Neverthe- 
less, one can still see that a naive perturbative expansion 
in B does not always work. In the naive expansion, one 
finds 



AE^^—B 
2 



(12) 



with given by Eq. ([5]), which recovers the conven- 
tional result. However, around the (in)stability threshold 
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Figure 3: The real part of the energy shift AE as function 
of for a fixed magnetic field strength B. The red (solid) 
curve is obtained from Eq. (|9]) while the blue (dotted) from 
Eq. (|12|) and Eq. ([5|. The parameters are chosen r — 0.9 and 
B = 0.05. 



Afn. = M + m, the naive expansion breaks down, as can 
be seen from Fig.[3]where we plot the energy shift Eq. ^ 
compared to the result of the naive perturbative expan- 
sion: Eq. with Eq. (O . It is clear that the two results 
are very different around the threshold which here is at 
/i = 0.1. The size of the region where the two results are 
different is proportional to the strength of the magnetic 
field. 

In Fig. 0] we again compare the perturbative and 
non-perturbative results, but now as a function of the 
magnetic-field strength. The masses are fixed such that 
the particle is stable for solid and long-dashed curves 
and unstable for medium- and short-dashed curves. In 
either situation there is a kink appearing at some value 
of the magnetic field, which indicates the crossing over 
the decay threshold. When is quasi-stable, /i + r = 1, 
the kink appears at i? = 0, which makes it impossible to 
define the moments as derivatives of the energy response 
with respect to the e.m. field. 

Integration over the Feynnian-parameter x in Eq. ^ 
yields more insight into the non-analytic dependence on 
the e.m. field. The result can be written as 



AE 



(13) 

where i? is non-analytic in B : 
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A In 



(a + A)(/3 + A) 
(a-A)(/3-A)' 
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Figure 4: The real part of the energy shift AiJ as function of 
the magnetic field B for a fixed [i. The parameters are chosen 
as r = 0.9. The red (solid) curve is obtained from Eq. (j9l) 
and the blue (long-dashed) from Eq. (fT^ with Eq. (O for 
^i, — 0.09. The brown (medium-dashed) curve is obtained from 
Eq. ((9)1 while the purple (short-dashed) curve from Eq. (|f 2|) 
for = 0.11. 



with 



1 



a 



1 + - - B] , 



2{1-B) 

P = rrr^ {l-r^ + ^l^-B] , 



1/2 



2(1-5) 

A = p-rV(l-S) 

while the analytic terms are contained in 

A = -2 + /31n/i^ + alnr^ 
/i2(l - In ^2-) _ ^2(^ 



(15) 



In r^) 



2{a + r){l-B) 



(16) 



(14) 



From the expression for we can readily see that a Tay- 
lor expansion in B only make sense when the condition 
of Eq. ([3]) is satisfied. 

The masses of particles are rarely tuned to the extent 
that the condition Eq. ^ is in danger. One field of 
applications where one does need to pay attention is lat- 
tice QCD. In modern lattice studies the e.m. moments of 
hadrons can directly be accessed using the background 
e.m. field method [8|. However, the field strength can- 
not be arbitrarily small, the periodicity condition poses 
a lower bound. In the case of magnetic field the bound 
is: eB > 2TT/{a'^L), or in best case il: eB > 27r/(aL)2, 
with length a and integer L being respectively the lattice 
spacing and size. For typical modern lattices the lowest 
possible value of the magnetic field can be as large as 10^'* 
Tesla. Certainly in such strong e.m. fields the problem 
raised here becomes relevant and should be studied on a 
case-by-case basis. 
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Figure 5: The condition Eq. ([3]) for the A-nucleon-pion 
system, \ 2M^{M^-Mt,-m^) \ < 1' plotted for the range of 
fields used in \eB\ = 0.00108/a^ . . . 0.00864/a^ with 

1/a = 2 GeV as function of the pion mass. Red (solid) curve 
corresponds to the stronger field and the blue (dashed) to the 
weaker field. The Delta-nucleon mass difference is taken to 
be constant: Ma - Mjv = 0.293 GeV. 



the Standard Model (prior to the top-quark discovery) 
or some of the hadrons in chiral effective-field theory, 
reflect only the limitation of the calculational tech- 
nique. When the mass of the particle is near a decay 
threshold (quasi-stable state), a small external e.m. field 
may induce the decay instead of interacting with the 
particle's e.m. moments. We have formulated an exact 
condition for this effect to occur. In this situation an 
extra care should be taken in defining and determining 
the moments, as has been described in this work. The 
present and future lattice QCD calculations of hadron 
e.m. moment using the background e.m. field technique 
are a very likely subject to this problem. 
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One typical example would be the case of the A(1232) 
isobar, which magnetic moment has recently been com- 
puted u sing th e background field method for various pion 
Figure [5] shows how the condition 



masses 



Hi 



eB 



2Ma{Ma - Mn - mn) 



< 1 



(17) 



can be violated in this type of studies, but of course 
for very specific values of pion mass and the background 
magnetic field. We emphasize that the actual parameters 
in llOL do not violate the above condition, mainly 
thanks to the large values of pion mass used in these 
works. However, current lattice calculations begin to ap- 
proach the pion-mass range where this condition would 
be violated. It would be interesting to see how the non- 
analytic B^/^ behavior emerges in these calculations. Of 
course one can expect this behavior to be shielded by the 
finite volume effects, the question is to which extent. 

To conclude, the singularities found in calculations 
of the e.m. moments of particles, such as W-boson in 
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